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Since 2008, many new families of iron-based high tem- 
perature (high-r c ) superconductors have been discovered 
QH4]. Unlike all parent compounds of cuprates that share 
a common antiferromagnetically (AF) ordered ground 
state, those of iron-based superconductors exhibit many 
different AF ordered ground states, including collinear-AF 
(CAF) state in ferropnictides |5|, bicollinear-AF (BCAF) 
state in 11-ferrochalcogenide FeTe (5J [7), and block-AF 
(BAF) state in 122-ferrochalcogenide £0.8^1.6^2 (H- 
While the universal presence of antiferromagnetism 
suggests that superconductivity is strongly interrelated 
with magnetism, the diversity of the AF ordered states 
obscures their interplay. Here we show that all magnetic 
phases can be unified within an effective magnetic model. 
This model captures three incommensurate magnetic 
phases, two of which have been observed experimentally. 
The model characterizes the nature of phase transitions 
between the different magnetic phases and explains a 
variety of magnetic properties, such as spin-wave spectra 
and electronic nematism. Most importantly, by unifying 
the understanding of magnetism, we cast new insight on 
the key ingredients of magnetic interactions which are 
critical to the occurrence of superconductivity. 

The iron-based superconductors can be divided into two 
major classes, ferropnictides and ferrochalcogenides. Su- 
perconductivity in both classes, like cuprates, arises from 
electron- or hole-doping of their AF parent compounds. 
Hence, determining magnetic interactions in the parent com- 
pounds is extremely important in the identification of key 
magnetic interactions relevant to the development of super- 
conductivity. 

However, because of the diversity of the magnetic or- 
ders and the fact that the magnetic properties exhibit the di- 
chotomic behavior of both local moment and itinerant elec- 
tron aspects, it has been extremely difficult to find a consistent 
magnetic model to describe the magnetism of iron-based su- 
perconductors. Theoretically, magnetism can be explained by 
either local moment models where local spins interact with 
each other or itinerant electron models where magnetic or- 
der arises from nested Fermi surfaces. The former is appro- 
priate for insulating materials such as the parent compounds 
of cuprates while the latter is suitable for metallic systems, 
such as chromium. However, the iron-based superconductors 
include rather diversified materials whose parent compounds 
can have either metallic or insulating ground states. More 
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specifically, the parent compounds of ferropnictides are bad 
metal 1 1 ] while the newly discovered 122-ferrochalcogenide, 
^0.8^1.6^^2 is a block AF insulator (8). Moreover, it was 
shown that even in the insulating parent compounds, itinerant 
electron aspects have to be included [10 ] and in the metallic 
parent compounds, local moment aspects are manifested I'll]. 
Such a dichotomy of the magnetic properties lfT0l - FT6Tl leads to 
many diversified viewpoints on what is the proper model to 
describe the magnetism fT7U22ll . 

In ferropnictides, the CAF state as shown in fig 1(b), in prin- 
ciple, can be explained by Fermi surface nesting between hole 
pockets at T and electron pockets at X l23l because they are 
connected by the CAF ordered wavevector Qcaf = (0, n) in 
unfolded reciprocal space. However, the mechanism fails to 
explain the BCAF and BAF states shown in figl(c,d). In these 
two cases, no matching Fermi pockets can be connected with 
their ordered wavevectors Qbcaf = (Vr/2,7r/2) and Qbaf = 
On 1 5,7i 1 5). Within a local moment picture, the CAF state 
can be naturally obtained in a Heisenberg model with near- 
est neighbor (NN) J\ and next nearest neighbor (NNN) J2 
magnetic exchange interactions if J 2 > |/il/2 lH71fT9l . How- 
ever, this model is quantitatively incompatible with the large 
anisotropy of the NN exchange interactions along the ferro- 
magnetic (FM) and AF directions in the CAF state determined 
in neutron scattering experiments ifTTIl . Intriguingly, recent 
neutron experiments in ferrochalcogenides show that the sim- 
ilar large NN anisotropy exists in both BAF lfT2l and BCAF 
states |[T0l . Moreover, there is almost no anisotropy between 
the NNN exchange interactions even though the magnetic 
configurations along the two NNN coupling directions in both 
states are different (l0l[T2l as shown in figl(c,d). Two com- 
pletely different solutions backed by electronic structure cal- 
culations have been proposed to solve the incompatibility in 
ferropnictides: one emphasizes orbital order l24l]26l and the 
other suggests biquadratic spin interaction terms l20l . How- 
ever, both solutions are inadequate in understanding the BAF 
and BCAF phases. For example, electronic structure calcula- 
tions and orbital ordering mechanism also suggested an large 
NNN anisotropy in ferropnictides [ 27, 28 ], which is inconsis- 
tent with experimental results fT0l[T2ll . 

Here we attempt to formulate a minimum effective spin 
model that unites the description of the magnetic properties 
of the parent compounds of the different classes of iron-based 
superconductors. The model has to preserve the tetragonal 
lattice symmetry so that it is capable of providing us the 
detailed relations between different magnetically ordered 
states as consequences of spontaneous symmetry breaking at 
low temperature. The model should be able to capture all the 
magnetically ordered ground states observed in iron-based 
superconductors, to explain their correct spin- wave spectra 
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and the anisotropy of magnetic exchange interactions, and 
to predict possible new states including incommensurate 
magnetic states. In the following, we will show by including 
the NN biquadratic interaction term as proposed in ref. l20lL 
but not the NNN biquadratic interaction term, and the next 
next nearest neighbor (NNNN) AF Heisenberg interactions 
J 3 EH 130) in J x - J 2 - J c model EU, we can fulfill above 
requirements. 

Model Hamiltonian We start with the following gen- 
eral Hamiltonian, 

H — ^ ^ [Jjj$ i^ n • S — Kij(S • S j^ n ) ] + J c S i iH • S i, n +l (1) 
ij,n i,n 

where J t j describe in-plane magnetic exchange interactions, 
J c is inter-plane magnetic coupling along c-axis (between iron 
layers) and K t j are in-plane non-Heisenberg biquadratic cou- 
plings. In the minimum model proposed here, we choose non- 
vanishing Jij = J i , J2 or J3 if and only if i, j are two NN, NNN, 
or NNNN sites respectively and K t j = K if and only if i, j are 
two NN sites. The interactions are sketched in fig 1(a) by the 
dashed lines. 

We note that the model is a natural extension of models 
in ref.OjO |20l |29l [30) proposed for ferropnictides and fer- 
rochalcogenides before. However, all previous models only 
describe particular family and fail to provide a comprehensive 
understanding of different magnetic states. 

Exact Classical Phase Diagram The classical phase 
diagram of the model can be obtained exactly. In fig2, we 
draw a typical phase diagram in the J3/\J\\-J2/\Ji\ plane by 
taking KS 2 /\J\\ = 0.2. The phase diagram is almost symmet- 
ric between J\ > (the right part of fig2) and J\ < (the left 
part of fig2). For J\ > 0, there are three commensurate phases 
labeled as AFM, CAF and BCAF in fig2, which exactly 
describe the static magnetic states of the parent compounds 
of curpates, ferropnictides and 11-ferrochalcogenide FeTe 
respectively. There are also two incommensurate phases 
sandwiched between the commensurate phases with ordered 
incommensurate wavevectors (q,n) or (7r, q) (labeled as IC1) 
and (q, q) (labeled as IC3) respectively. The static (q, q) phase 
was observed in Fe\ +y Te when y > 0.1 [6]. Although no 
static (q,n) phase has been detected, the (q,n) incommen- 
surate spin fluctuations have been observed in FeTe\- x S e x 
OTI [32l , electron-overdoped Ba{Fe\- x Co x \As2 1 33 ], and 
hole-doped KFe^As-i l34l . If J\ is switched to negative, 
namely ferromagnetic (FM), the AFM phase becomes a FM 
phase and (q, n) becomes (q, 0) or (0, q) (labeled as IC2). 

More specifically, we can determine phase transition 
boundary. We scale other parameters with J\ as K = KS 2 /J\, 
3 2 = JilJ\, and ^3 = /3//1 for simplicity. The phase 
boundary between the BCAF and the (q, n) incommensurate 
phase is determined by 4(J 3 - \) 2 - (J 2 - \) 2 = (K - \ f 
which defines the upper branch of a hyperbolic curve centered 
at 0/2,^3) = (1/2,1/4). The phase boundary between the 
AFM and BCAF phases is determined by J3 = -\Ji + \ 
and that betwee BCAFM and CAFM phases is determined 
by 1/3 = \J2. The incommensurate states appear only when 
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TABLE I: The values of the magnetic exchange interactions in J\ - 
J 2 -J 3 -K model obtained from the experimental results of different 
parent compounds of iron-based superconductors 1 10 -12 16 1. 



K < 0.5. We emphasize that finite positive K > (1.5 - V2) 
and 1/3 > 0.25 are necessary conditions for the appearance 
of the BCAF phase. The incommensurate wavevectors 
of the three incommensurate phases can also be explicitly 

2/2+1 . 



determined: (1) (g, 0) or (0, g) phase with q = arccos 
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(2) (q,n) or (7r, q) with q = arccos ^-Jp; and (3) (q, q) with 
q - arccos — J— - , where A - K - 2/3 . 

2{J2-A) 

Commensurate Phases The model captures three commen- 
surate phases: AFM, CAF and BCAF. In these commensurate 
phases, the biquadratic interaction term effectively creates 
the anisotropy of the NN magnetic exchange interactions 
by taking a meanfield decoupling EDI . Depending on the 
spin alignment of two NN sites, J\ a = J\ + 2KS 2 if it is AF 
and J\t = J\ - 2KS 2 if it is FM. Therefore, effectively, our 
model becomes a J\ a - J\b - Ji - J3 - J c model in these 
phases. Experimentally, the spin wave excitations in the 
parent compounds, CaFe2As2 ifTTTl and BaFe2As2 fT6lL were 
fitted well to the J\ a - - J2 - J c model, and those of FeTe 
was fitted well to a J\ a - J\b - J2 - J3 - J c model 1 12 ]. 

More interestingly, in the ferrochalcogenide KQ ^Fe\^Se2, 
featuring intrinsic vacancy ordering, the fitting of spin wave 
excitations to a J\ a - - J2 - J' 2 - J3 model has concluded a 
large FM J\ a , AF and very small anisotropy between two 
AF NNN couplings J2 an J' 2 \ 10]. Our model hence success- 
fully describes the magnetism of this new material as well. 
The presence of vacancy ordering does not drastically affect 
the magnetic interactions. Instead, it reduces the magnetic 
frustration to stabilize the BAF order f[35l . 

From the experimental results of spin wave excitations |[T0l - 
[l2l[T6l, we now can extract the magnetic exchange parameters 
of our model for different parent compounds. The results are 
summarized in Table 1. We note that FeTe is near the bound- 
ary of the BCAF phase and incommensurate phases. The val- 
ues listed in Table 1 is within the error bar of experimental 
values in [12]. This table displays a central message that all 
iron based superconductors share a similar AF NNN exchange 
interaction J2. However the sign of J\ is different between 
ferropnictides and ferrochalcogenides and a significant AF J3 
exists in ferrochalcogenides but not in ferropnictides. 



Spin Excitations in Incommensurate Phases The model 
also predicts three incommensurate phases. We can calculate 
the distinct features of spin excitations in these three phases. 
The typical spin excitation spectra in the three incommen- 
surate phases as the function of wavevectors S are shown in 
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fig3(a)-(f). In general, spin excitations include two branches, 
which can be identified as an acoustic mode and an optical 
mode. In the (n, q) phase, the two modes connect with each 
other at two new incommensurate positions 8 = ±qo when the 
wavevector (n, 8) varies as shown fig3(a). If we calculate the 
intensity of these spin wave excitations measured by neutron 
scattering, an hour-glass like behavior along this direction 
becomes prominent at low energy as shown in fig4. This 
behavior has been recently reported in FeS ^0^^0^ 1361 and 
a clear explanation was not given before. The dispersion 
along (8,n) as varying 8 has much larger energy dispersion 
than the one along in, 8) and reaches maximum at (n,7i) 
which is consistent with experimental results observed in 
ref. (31) as shown in fig3(b). In (0, g) phase, the dispersion 
of spin excitations along (0, 8) as varying 8 is very similar to 
the one along (8,n) in (n,q) phase as shown in fig3(c) and 
that along (8, 0) is similar to the one along (n, 8) as shown 
in fig3(d). In the (q,q) phase, the two spin wave modes 
connect at (0,7r) point as shown in fig3(e,f) and fig3(g,h), 
displaying a Dirac-type dispersion. The two modes has much 
larger energy separation in the case of J\ > (shown in 



fig3(e,f) ) than J\ < (shown in fig3(g,h). More detailed 
spin wave properties are included in supplement materials. 
These distinct features can be used to determine the effective 
magnetic exchange couplings even if the incommensurate 
order is not static. 

Nematism and Effective Field Theory Both the CAF 
and BCAF states break the C 4 rotational symmetry of the 
tetragonal lattice. The rotational symmetry breaking can be 
separately described by an Ising or nematic order as shown in 
ref. EH [19). Without the specific biquadratic term K, when 
the parameters of the J\ - Ji - J c model are fixed in the 
CAF phase region, a weak biquadratic term can be developed 
through the 'order by disorder' mechanism |Qj3 UU [33 and 
the nematic phase transition can take place at a transition 
temperature T N higher than the CAF transition temperature 
T c if the inter-layer coupling J c is much weaker than J 2 lfT8) . 
This physics can be analytically described in the continuum 
limit. As show in ref. [18), the effective field theory of the 
J\ - Ji - J c model in the continuum limt is given by 



CAF 



f d 2 Y ^ ^/2|V0 w , a (r)| 2 - J c $ n ,a( V ) ' </>n + l,a( V )\ -g^j ( r ) ' feWf + Jl X $n,l(X)d X dy4>n,2(X), (2) 



r 



where we use the same notions as ref. [18]: $ n ,a=i,2 spec- 
ify the two AF Neel orders in the two sublattices of the 
tetragonal lattice shown in fig 1(a) ( for simplicity, we take 
S = 1 in this section). The nematic order is defined to 
be c — 2g(0 w j(r) • (p n ^))- Without the biquadratic term 
K, g ~ 0.13/j7/2- With this term, we just need to modify 
g ~ 0.13/j7/2 + K. Therefore, the calculations and the results 
in ref. |[T8l are still valid. The large g value due to the spe- 



cific biquadratic K term simply enhances the nematic order 
and increases the temperature range of spin nematic fluctua- 
tion above T N , which has been observed experimentally lfT6l . 

In the case of BCAF states, the similar effective model for 
J\-Ji — h model has also been derived in ref. (38) • Including 
J c , the effective field theory of the J1-J2-J3- J c model close 
to the BCAF states can be written as 



Hbcaf = f (i 2 r^|i/ 3 |V^(r)| 2 -/ c ^(r)-0 w+1?( ,(r) - g' ^{[<^,i(r) • <^, 3 (r)] 2 + [$ n ,i(x) ' 0M( r )] 

n,a n 

-Ji J d 2 r ^ ^ n ,i(r) • V x w , 2 (r) + $ n ,4(r) ' V x w , 3 (r) - 0„, 2 (r) • V^ n>3 (r) - w ,i(r) • V y $ nA (r) 



(3) 



r 



where <p nfil with a - 1, 2, 3, 4 are four Neel order parameters 
defined in the four sublattices of the tretragonal lattice spec- 
ified by the J3 coupling and g' ~ 0.1 3/|//3. Differing from 



the case of the CAF, the biquadratic term in the case of the 
BCAF adds an additional coupling term into above effective 
Hamiltonian, which is given by 



H K = -K J d 2 rYj{$n,i(x) • fe(r)] 2 + [&,i(r) • $ nA (r)f + [^(r) • ^(r)] 2 + [^(r) • M (r)] 2 } (4) 
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The complete field description of the J\ - J2 - J3 - J c - K 
model near the BCAF phase is given by H T = H B caf + H K . 
In this model, there are four independent Ising or- 
ders (Ti = g'$ n ,\(X) ' 0«,3( r )' °~2 = g'4>n,l(X) ' 0«,4( r )> 

cr+ = K($ nt i(r) + 0„, 3 (r)) • (0„, 2 (r) + <^, 4 (r)) and 
cr- = K($ nA (r) - w , 3 (r)) • (0„, 2 (r) + 0„, 4 (r)). A detailed study 
of this effective field model will be present elsewhere. Here, 
we simply discuss the case when K is much larger than g' 
to demonstrate that the effective field model indeed captures 
phase transitions. In this case, the transition between the 
BCAF state and the (q,q) incommensurate state is controlled 
by the Ising order <x+. It is easy to show that when |<x+| is 
large, the ground state is the BCAF state. The transition takes 
j 2 

place at |<x c | = When \cr+\ < |<x c |, the model is in an 
incommensurate (q, q) state with q = j^(\o~ c \ - \o~+\). 

Discussion By describing the magnetism of the differ- 
ent parent compounds of iron based superconductors in a 
single effective magnetic model, we can cast new insight on 
the microscopic origin of magnetism. From the magnetic 
exchange coupling paramters of the effective model, it is 
very clear that the magnetism is neither purely local nor 
purely itinerant, rather it is a complicated mix of the two. 
The presence of significant NNNN coupling J 3 suggests 
the local exchange mechanisms such as superexchange or 
double exchange are not enough to account for all magnetic 
exchange interactions. Moreover, the sign change of J\ 
between ferropnictides and ferrochacogenides suggests 
that the NN exchange interactions are sensitive to subtle 
difference in band structures. However, the robustness of 
NNN J2 interactions indicates that the NNN J2 coupling is 
most likely determined by local superexchange mechanism. 
In ref. ||20L the authors suggest that the NN biquadratic term 
K in ferropnictides stems from the strong magnetostructural 
coupling. However, this explanation does not provide an 
understanding of the absence of NNN biquadratic term in the 
BCAF state of FeTe since the lattice distortion in FeTe is 
monoclinic rather than orthorhombic as in ferropnictides! 6]. 
More throughout studies of the origin of the biquadratic term 
are needed. 

The model reveals the significant difference between fer- 
ropnictides and ferrochacogenides: the sign difference of J\ 
and the large AF J 3 in ferrochcogenides. These significant 
differences may suggest the importance of the p orbitals of As 



or Te/Se on the influence of magnetism. So far, most theo- 
retical models are constructed based on the d orbitals of irons 
with onsite interactions. Since the effect of electron-electron 
correlations is believed to be weaker in ferropnictides than 
in ferrochalcogenides, one would expect the range of mag- 
netic interactions should be shorter in ferrochalcogenides than 
in ferropnictides, which contradicts the existence of large J3 
interactions in ferrochalcogenides but not in ferropnictides. 
This contradictory can be resolved if the significant parts of 
magnetic exchange interactions are generated through the p- 
orbitals of As or Se/Te. The effective magnetic exchange in- 
teractions obtained from onsite electron-electron interactions 
are not enough to account for entire magnetic exchange cou- 
plings. This also explains that why magnetism is so sensitive 
to the distance of As or Se/Te away from iron planes l39l be- 
cause the distance may strongly affect the mixture of p orbitals 
in electronic structure. 

The model, as an effective low energy model of magnetism, 
also tells us the power and limitation of LDA calculations 
performed for iron-based superconductors where the effect of 
electron-electron correlation can not be ignored. Without any 
doubt, the LDA calculations explain many magnetic proper- 
ties in iron-based superconductors. For ferropnictides, the 
LDA results of J\ a , J\b and J 2 values are in a good agreement 
with experiments l22l . However, LDA calculation wrongly 
predicted the large anisotropy of J 2 in the BCAF and BAF 
states (ST] [28). This failure is not surprising since the LDA in 
magnetically ordered state is simply a complicated meanfield 
approach. 

It has been shown that the high energy magnetic excita- 
tions in electron-doped ferropnictides are very similar to those 
of parent compounds |40l . This proves that the short range 
magnetic correlations in superconducting states are still dom- 
inated by the magnetic exchange interactions determined in 
the corresponding parent compounds. The doping destroys 
the long range magnetic correlation but not the short range 
interactions. Especially, the J 2 magnetic exchange interac- 
tions should be expected to vary little against doping. In 
FeTe\- x Se x , the incommensurate spin excitations are rather 
robust against the replacement of Te by SeEU [32|. This fact 
suggests that J3 is relatively stable against the replacement in 
this family of materials as well. Therefore, if AF exchange 
couplings are responsible for superconductivity, we expect 
both J 2 and J3 plays a significant role in superconductivity 
of ferrochacolgenides. 
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FIG. 1: (a). The sketch of magnetic interaction parameters of J1—J2-J3- K. (b) The collinear-antiferromagnetic state (CAF) in ironpnictides, 
for example CaFe 2 As 2 . (c) The bicollinear-antiferromagnetic state (BCAF) in FeTe. (d) The block- antiferromagnetic state (BAF) with 
V5 x V5 vacancy ordering in K ^Fe^ 6 Se 2 . 
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FIG. 2: The classical phase diagram of the J\ - J 2 - J 3 - K model at K = 0.2 J \ . There are total four commensurate and three incommensurate 
magnetic phases. Their labels and the ordered wavevectors are specified in the figure. 
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FIG. 3: The spin waves in incommensurate phases: (a) along (l,6)n and (b) along (6, l)n in the (tt, q) phase, (c) along (0,6)7T and (d) 
along (6, 0)tt in the (0, q) phase, and (e)/(g) along (8, 8)tt and (f)/(h) along (5, -5)tt in the (q, g)-phase. The parameters are fitted to be 
S = I, (J u J 2 /\JilJ 3 /\JilKS 2 /\J 1 \) = (1,0.6,0.06,0) in (a) and (b), (-1,0.6,0.06,0) in (c) and (d), (-1,0.6,0.5,0.05) in (e) and (f), and 
(1,0.6,0.5,0.05) in (g) and (h). 
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